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We use methods of quantum field theory in toroidal topologies to study the A'^-component D-dimensional 
massive Gross-Neveu model, at zero and finite temperature, with compactified spatial coordinates. We discuss 
the behavior of the large- A'^ coupling constant (g), investigating its dependence on the compactification length 
(L) and the temperature (T). For all values of the fixed coupling constant (A), we find an asymptotic-freedom 
type of behavior, with g —>■ as L —>■ and/or T — >■ oo. At T = 0, and for A > Ac (the strong coupling 
regime), we show that, starting in the region of asymptotic freedom and increasing L, a divergence of g appears 
at a finite value of L, signaling the existence of a phase transition with the system getting spatially confined. Such 
a spatial confinement is destroyed by raising the temperature. The confining length, Lc , and the deconfining 
temperature, TJ , are determined as functions of A and the mass (m) of the fermions, in the case of Z3 — 2,3,4. 
Taking m as the constituent quark mass (~ 350 MeV), the results obtained are of the same order of magnitude 
as the diameter (~ 1.7 fm) and the estimated deconfining temperature (~ 200 AdeV) of hadrons. 

PACS numbers: ll.lO.Kk; ll.lO.Wx 
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I. INTRODUCTION 

The strong interaction among quarks and gluons, the con- 
stituents of the hadronic matter, has such a structure that ob- 
ligates them to live spatially confined, at low temperatures, 
within distances ^ 1 fni in colorless states. It is usually ac- 
cepted that in an earlier stage of the Universe, as it cooled 
down, quarks and gluons condensed into hadrons at an esti- 
mated temperature of the order of 200 MeV. At very high en- 
ergies, deep inelastic scattering indicates that the quarks are 
nearly free, a regime denominated as asymptotic freedom. 

In the standard model, the theory of strong interactions is 
quantum chomodynamics (QCD), which should describe such 
facts, accounting also for the nuclear forces. However, QCD 
has a very involved mathematical structure, which practically 
prevents us from finding analytical results taking into account 
both confinement and asymptotic freedom. Lattice calcula- 
tions have been implemented to simulate the behavior of the 
theory in the confining region, both at zero and finite tem- 
perature, providing (among other results) an estimate of the 
deconfining temperature. Rigorous QCD calculations, both at 
zero and finite temperature, have been worked out [1,2] but 
mainly treating the asymptotically free domain at high ener- 
gies or high temperatures, where perturbation theory is appli- 
cable. 

Due to the difficulty of treating QCD analytically, phe- 
nomenological approaches and studies of effective, simplified, 
models have been stimulated along the years to give clues 
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to the behavior of hadronic systems. A celebrated effective 
model, which shares with QCD some basic properties, is the 
Nambu-Jona-Lasinio (NJL) model yfl. One of its sectors, pro- 
viding the simplest effective model which may be considered 
as describing quark interactions, is a direct four-fermion cou- 
pling, where gluon fields and color degrees of freedom are in- 
tegrated out, resembling the Fermi treatment of the weak inter- 
action. This corresponds to the Gross-Neveu (GN) model |3l, 
considered in space-time dimension D = 4. Although the 
GN model is not renormalizable for dimensions greater than 
D = 2, the Euclidian model has been shown to exist and has 
been constructed for I? = 3 in the large- A^ limit p]. But, 
within the spirit of effective theories, perturbative renormal- 
izabilityis not a requirement to have a physically meaningful 
model IMS- 

The GN model, as a prototype model for intera cting 
fermions, has been analyzed extensively in recent years lllll - 
V(J], including the study of continuous and discrete chiral sym- 
metry 1I2U I22I1 . For the version with A^ massless fermions 
in (2H-l)-dimensions, for instance, a chiral symmetry break- 
ing is found in perturbative analysis, with the restoration 
of such a symmetry at finite temperature (J231. In particle 
physics, these results provide insights into the intricate struc- 
ture of the hadronic matter, such as for the quark confine- 
ment/deconfinement phase transition 11241 [2511 . 

The four-point contact interaction of the GN model is simi- 
lar to the delta interaction in the BCS theory of superconduc- 
tivity. In the latter case, as in other systems of condensed mat- 
ter, the susceptibility arising from the linear response theory 
has a divergence at a finite temperature indicating the exis- 
tence of a second-order phase transition between a disordered 
and a condensed phase [26] . Spontaneous symmetry breaking 
is the common feature underlying all these phenomena |27]. 
Recently ll28ll . we have shown that such instability appears in 
the one-component massive tridimensional GN model at fi- 



nite temperature. Here, we intend to use similar treatment to 
investigate the existence of a phase transition in the massive, 
iV-component, GN model. 

We employ methods of quantum field theory in toroidal 
topologies ll29l - [34ll to extend previous results [i35., ,36.1 for the 
Euclidian massive, A^-component, GN model in the large- 
N limit. This amounts to consider the GN model in a D- 
dimensional space-time and to compactify a rf-dimensional 
(d < D) subspace. The compactification is a generalization 
of the Matsubara procedure. The Matsubara imaginary-time 
formalism corresponds to considering fields in a space with 
topology S^ X M-'-'^^, where §^ is a circumference of length /3, 
with periodic (anti-periodic) boundary conditions for bosons 
(fermions). Such a compactification of the Euclidian time 
can be directly generalized to include the compactification of 
space coordinates as well. This allows us to consider field the- 
oretical models with spatial constraints, at zero or finite tem- 
perature, by using generating functionals with a path-integral 
foi-malism on the topology S^i x • • • x S^" x M^"'^ L29.-.32J. 
These ideas have been established recently on a firm founda- 
tion [3J, 3?) and applied in different physical situations, for 
example: for spontaneous symmetry breaking in the compact- 
ified 0^ model Il37l - [39ll : for second-order phase transitions 
in superconducting films, wires and grains rw-421; for the 
Casimir effect for bosons and fermions [43-48|; for size ef- 
fects in the NJL mo del ll49l - l53[l ; and, for electrodynamics with 
an extra dimension \5^ . 

We treat particularly the cases D = 2,3,4 with all spa- 
tial dimensions compactified, initially at zero temperature and 
then discuss temperature effects by compactifying the imagi- 
nary time in a length (3 = T^^, T being the temperature. This 
corresponds to considering the system contained in a paral- 
lelepiped box, with anti-periodic boundary conditions on its 
faces, at finite T. We study the behavior of the system as a 
function of its size and of the temperature, concentrating on 
the dependence of the large- A^ coupling constant on the com- 
pactification length and temperature. We show that, even at 
T = 0, a singularity in the 4-point function may appear driven 
by changes in the compactification length suggesting the ex- 
istence of a second-order phase transition in the system. This 
can be interpreted as a spatial confinement transition, which 
may be present in the massive version of the GN model. 

We use concurrently dimensional and analytic regulariza- 
tions and employ a subtraction scheme where the polar terms, 
arising from Epstein-Hurwitz generalized zeto-functions, are 
suppressed. Results obtained with this procedure have similar 
structure for all values of D, which gives us confidence that 
they are meaningful for the 4-dimensional space-time. This is 
reinforced a posteriori by the fact that the numerical results 
found for Z? = 4 are of the same order of magnitude as the 
corresponding values for D — 2 and D = 3. For the massive 
GN model in the large- A^ limit, discussed in the present paper, 
we obtain simultaneously asymptotic freedom type of behav- 
ior and spatial confinement, in the strong coupling regime, for 
low temperatures. We also show that, as the temperature is 
increased, a deconfining transition occurs. We calculate the 
values of the confining lengths and the deconfining tempera- 
tures for D = 2, 3, 4. 



This article starts by discussing, in Sec. 2, the D- 
dimensional massive, A^-component, Gross-Neveu model 
with d (< D) compactified dimensions. In Sec. 3, the calcu- 
lation of the effective large- A^ coupling constant is carried out 
for the cases with D = 2, 3, 4 at zero temperature. Tempera- 
ture effects are presented in Sec. 4. The last section provides a 
comparison of our estimated confining lengths and deconfin- 
ing temperatures with experimental values. 



II. COMPACTIFIED GROSS-NEVEU MODEL 

The massive GN model in a _D-dimensional Euclidean 
space is described by the Wick-ordered Lagrangian density 

— u — 

C =: ib{x){if + m)i^{x) : +-(: ^{x)^{x) -.f , (1) 

where m is the mass, u is the coupling constant, x is a point 
of M^ and the 7's are the Dirac matrices. We consider the 
GN model in its A^-component version, so that "ipi^) repre- 
sents a spin i field having A^ (flavor) components, '4>"-{x), 
a = 1,2, ...,N, with summations over flavor and spin indices 
being understood in Eq. ([T]l. We calculate quantities of inter- 
est by taking the large- A^ limit where N —> 00 and w — >■ in 
such way that Nu — A remains finite. Throughout the text we 
use natural units with h ~ c = fcs = 1. 

Our main goal is to determine the large- A^ (effective) cou- 
pling constant when d (< D) Euclidian coordinates, say 
xi,...,Xd, are compactified, that is, considering the sys- 
tem in a topology S^^ x • • • S^'' x K^"''. This corresponds 
to restricting the coordinates Xi to segments of length Li 
(i — 1, 2, ....d), with the fermionic field iIj{x) satisfying anti- 
periodic boundary conditions. If all Xi are spatial coordi- 
nates, the model refers to the system compactified in a d- 
dimensional box at zero temperature while, with one coor- 
dinate being the Euclidian time (say Xd), one has the system 
with d — I compactified spatial dimensions at finite temper- 
ature; in this latter case, Ld would stand for (3 — 1/T, the 
inverse of the temperature. For massless fermions, this spa- 
tial compactification with anti-periodic boundary conditions 
is equivalent to considering the system constrained to "live" 
inside a parallelepiped "box", with edges Li (i = 1,2, ....d), 
under bag model conditions (no outgoing currents) on paral- 
lel, opposite, faces [55., .56.1 . In our case, to calculate n-point 
functions, we apply the generalized Matsubara prescription, 
which amounts to modifying the Feynman rules performing 
the replacements 



h -^ V, = 



d^k 

(27r)^ 



F(k) 



27r(n, + i) 



Li . . .Ld 

+00 

X E 

{rii} — — 00 



i = 1,2, ...,d; 



(2a) 



jD-di 



(2^) 



(2b) 



where {rii} = {ni, . . . ,nd}, with n^ G Z, and k is a 
{D — d) -dimensional vector in momentum space; the dis- 
creet momenta Vi are referred to as Matsubara frequencies. 
Additionally, it should be pointed out that the choice of anti- 
periodic boundary conditions for the spatial compactification, 
instead of the simpler periodic ones, is also due to the fact 
that they emerge naturally in the generalization of the KMS 
(Kubo-Martin-Schwinger) conditions satisfied by correlation 
functions for fermionic fields in toroidal topologies 11331 13411 . 
We shall define the large- A^ effective coupling constant be- 



tween the fermions in terms of the 4-point function at zero 
external momenta. The {L^} -dependent four-point function, 
at leading order in -^, is given by the sum of chains of one- 
loop (bubble) diagrams, which can be formally expressed as 



rW(0;{L,},«) = 



l + NuIlDd{{L,})' 



(3) 



where the {Li} -dependent one-loop Feynman diagram is 
given by 



UDdiiU}) = 



Li ■ ■ ■ Ld 



E 



{ni} = -c 



jD^di 



(27r)^-'' 
I 



Ed 9 



k^ + Eti-? 



(4) 



Prior to defining an effective large-TV coupling constant, we the dimensionless quantities 6; = {mLi)~'^ {i = l,...,(i) 
have to deal with the ultraviolet divergences of I{od{{Li]). and qj — kj/2TTm {j — d + 1, . . . , D), in terms of which the 
To simplify the use of regularization techniques, we introduce one-loop diagram is written as 



J 



TlDdi{h})^TlDdis;{h})\^^ 



,D-2 



s=2 



^^^^Vbi---bd I ^Uodis; {h}) - Uodis - 1; {h}) 



s=2 



(5) 



where 



oo „ 

UDd{^^■,{b^})^ E j 



q' + E-=i^^-K- + 5)' + (2'^) 



(6) 



r 



We find from Eq. Q that liod has dimension of m^ ^, which transformed into sums over integers leading to 
is inverse of the mass dimension of the coupling constant. 



We employ a modified minimal subtraction scheme which 
uses concurrently dimensional and analytical regularizations. 
In this scheme, the subtracted terms are poles (for even D > 
2) of the Epstein-Hurwitz zeia— functions. First, using well- 
known dimensional regularization formulas to perform the in- 
tegral over q = {qd+i , • • ■ , Qd) in Eq. (|6]l, we obtain 



UDdi^J';{h}) = TT 



r(M) 
[zf(,7,6i,...,M 



4" 



E^^''(^' 



,...,46,,.. 



D-d I [Li — 



D-d \ 
1 ) 



r(M) 



d 

E^. 



E 



2 



(2^)- 



(7) 



The summations over half-integers in this expression can be 



+ Y. Zf (r;,...,4&„...,46„... 
-■■■+ i-lY Zf{^ Ah,,... Ahd) 



where K^ = n ^, 77 = u — 



D-d 



rj = 1^1 — "„ " and 



Zf{v,W^})-- 


00 

= >: 


d 




{n,} = -oo 


J='^ 



(8) 



(9) 



is the multiple (d-ditnensional) Epstein-Hurwitz zeta- analytically extended to the whole complex 77-plane 113911 . 
function. throu gh a g eneralization of the procedure presented in 



The Epstein-Hurwitz zeia-function Z'] (77, {a;}) can be Refs. ll57Ll58ll : then we find 



zf{rl,W^}) = 



■K2 



1 



d 00 



rU 






K„ 



2'Kh\ 



(10) 



where {cfq} represents the set of all combinations of the terms appearing in the parcels of Eq. (|8) and using the identity 
indices {l,2,...,(i} with 9 elements and Ka{z) is the 
Bessel function of the third kind. Consequently, the function 
Uodim {bi}) can also be analytically continued to the whole 
complex /i-plane. 

Taking Zj^ ivA'^i}) given by Eq. ( fTOl i. grouping similar 



N 

E 



-1 



iV! _ 1 
j\{N-j)\ "2^ 



(11) 



we obtain 

UDd(p;{b^}) = 

with WodilJ'; {bi}) given by 

WDd{p;{b^}) 



22m— 0,^2^- f I 



T{fi) V5i • • • bd 



2^-E2^^E E n 

J = l {P.}{cp,,=l,4} \fc=l 



r M- 



(-i)-p. 



^Pk 



D 



2-WDd{^i;{b^}) 



FDji^J.■,Cp,bp^,...,Cp^bp^), 



(12) 



(13) 



where {pj} stands for the set of all combinations of the indices {1,2, ... ,d} with j elements and the functions 

Fdj {jjl; ai , . . . , qj), for j = 1, . . . ,d, are defined by 




FDj{p;ai,...,aj) 



Substituting Eq. (fTSI) into Eq. ^ leads directly to an ana- 
lytic extension of ni3(i(,s; {bi}) for complex values of s, in the 
vicinity of s = 2. In fact, Hodis; {bi}) can be written as 

^Dd 'is) 



K„ 




(14) 



r 



nDd{s;{b^}) = np°'^^ 



(27r)^-2«+4r(s) 
x[2WDd{s■,{b^}) 
- {s-l)WDdis-l;{b^})], (15) 



where 



n 



polar 



(^) 



l^-^TT^is-l-D) 



r s- 1 



D 

'2 



(16) 



and the functions WodilJ'', {bi}) are given by Eq. ( fT3] l. We 
notice that the first term in this expression for IIud{s; {bi}), 



n^°^'*' (s), does not depend on parameters bi, that is, it is in- 
dependent of the compactification lengths Li (i = 1, . . . ,d}. 
At s = 2, because of the poles of the F-function, such a term 
is divergent for even dimensions D > 2. 



In order to obtain a finite single bubble function, we shall 
use a modified minimal subtraction scheme, where terms to be 
subtracted have poles appearing at the physical value .s = 2. 
Thus, the polar parcel given by Eq. (fT6t will be suppressed 
and, for the sake of uniformity, this term is also subtracted in 
the case of odd dimensions, where no pole of the F-function 
is present; in such a situation, we perform a finite subtraction. 
In this way, using the same notation, we define the finite one- 



loop diagram by the relation 



nn4{b^})- {Tin As; {b^}) -Tilths)} 



(17) 



Therefore, the finite one-loop diagram, which depends on 
the compactification lengths Li and arises from the regular 
part of the analytical extension of the Epstein-Hurwitz zeta- 
functions, is given by 



IlDdiM) = 



,D-2 



{2ny 



[2WDd{2;{h})-WDd{l;{b^})] 



(18) 



From now on, we shall deal only with finite quantities that 
are obtained following this subtraction prescription. Notice 
that, replacing bi by {mLi)^'^ in the above expression, we 
recover explicitly IV£)j^{{Li]). Now, we proceed to analyze 
the behavior of the large- A^ coupling constant in various cases. 



suggesting that the system presents an asymptotic-freedom 
type of behavior for short distances and/or for high temper- 
atures. 

From the extreme limits considered above two situations 
may emerge, as one changes the compactification lengths 
from to oo: either liod varies from oo to through pos- 
itive values, or liod reaches before tending to through 
negative values. The latter case, which may actually happen, 
would lead to an interesting situation where a divergence of 
the effective coupling constant would appear at finite values 
of the lengths Li. This possibility, and its consequences, will 
be investigated explicitly in the following subsections, con- 
sidering the compactified GN model at T = for space-time 
dimensions D = 2,3,4. The discussion of finite temperature 
effects is postponed to Sec.lIVI 



A. 2-D compactified GN model at T = 



III. COUPLING CONSTANT IN THE LARGE-TV LIMIT 

In field theories with four-fermion interactions, the cou- 
pling constant is defined in terms of the four-point function at 
fixed external momenta; here we choose p = 0. In this situa- 
tion, the coupling constant can be interpreted as measuring the 
strength of the interaction between the fermions. The large- A^ 
({&i} -dependent) coupling constant, for d (< D) compactified 
dimensions, is then obtained by substituting IiDd{{bi}) into 
Eq. ^ and taking the limit N -^ cx), m ^- 0, with Nu ~ A 
fixed; we get 



gDd{{bi\A) = jm NV^^\{Q,{hi\,u) 



Nu=\ 



A 



l + AHDdd&O)' 



(19) 



It is clear that, while gDd{{bi}^ A) depends on the value of the 
fixed coupling constant A in a direct way, its dependence on 
the compactifiation lengths is dictated by the behavior of liod 
as {bi] is varied. The dependence of god on {Li] and A is the 
main point to be discussed in the subsequent analysis. 

Limiting behaviors of the finite coupling constant god can 
be readily obtained from the fact that liod depends on {bi] 
through the Bessel functions of the third kind appearing in the 
definition of the functions Fjjj, Eq. ( fT4] i. First, if we let all 
the compactification lengths tend to infinity, that is {bi — > 0}, 
thus reducing the problem to the free space at T = 0, then 
n^d -^ and we obtain, consistently, that 



lim gDd[{bi],\) = A, 



(20) 



where A is the fixed coupling constant in free space at zero 
temperature. This is a consequence of the fact that Ki,{z — > 
oo) — > 0, for V integer or half-integer In the opposite limit, 
for any bi tending to cc (that is, if any compactification length 
Li goes to 0), the single bubble diagram liod — >■ oo, since 
A'i,(z) ^- as z ^- cx). This implies that the effective cou- 
pling constant gnd vanishes, irrespective of the value of A, 



For D ~ 2 and d — 1 (two-dimensional space-time with 
the spatial coordinate compactified), we put hi = {mL)~'^ in 
Eqs. (113114b . then Eq. ([T8]l becomes 



(21) 



n2i(i) = 2Ei{2jnL) - Ei{mL), 
where the function Ei (x) is defined by 

-. oo 

Ei{x) = - V [-Ko{xn) + {xn) Ki{xn)] . (22) 

n=l 

Notice that 11 is dimensionless for D = 2. 

We can calculate 1121 (i) numerically by truncating the se- 
ries appearing in the definition of the function Ei {y), Eq. ( |22] |. 
at some value n ^ M. For moderate and large values of niL 
(say, mL > 0.5), one can take for M a relatively small value; 
e.g., for mL — 0.5, choosing AI = 36 already leads to the 
correct value of 1121 to six decimal places. As mL increases, 
the value of M can be made smaller to give the same preci- 
sion. However, since the functions Ko{z) and Ki{z) diverge 
for z — ^ and the summation involves positive and negative 
parcels, the calculation of 1121 for small values of mi requires 
large values of A/; for m,L — 0.005, we need M « 4500 to 
obtain 1121 to six decimal places. Fortunately, the relevant be- 
havior of 1121 (L) appears for moderate values of m,L so that 
the numerical calculations are carried out in a short computa- 
tional time with very good precision. 

The function 1121 (i) is plotted as a function of mL in 
Fig. [T] From this figure and the numerical treatment of 
Eq. ( l2n i. we infer that 1121 (L) diverges (— > +oo) when L ^- 
and tends to 0, through negative values, as L — > oo. Also, we 
find that 1121 (L) vanishes for a specific value of L, which we 

(2) (2) 

denote by Li^[^, being negative for all L > L)^[^, and as- 
sumes a minimum (negative) value at a value of L denoted 



(2) 



for reasons that will be clarified later Numeri- 



by L 

cally, it is found that L 

and ng;'" 



(2) 



0.78 m-i,L 



(2) 



1.68 771 



-0.0445. This behavior of 1121 as L changes, 
particularly the fact that 1121 (L) < for L > l''^^^, leads 
to remarkable properties of the large- A^ coupling constant 



g2i{L, A). It is important to point out that such a dependence 
of the polarization on L is a direct consequence of the use 
of anti-periodic boundary conditions for the spatial compact- 
ification. Taking periodic boundary conditions (PBC), one 
would obtain n2i^'^(L) = Ei{rnL)/A which is positive for 
all values of L, and so no significant size-effect would exist. 



s 

0.2 




FIG. 2: Plots of the relative effective coupling constant, G — 
g2i{L, A) /A, as a function of mL for some values of A: 12.0 (dashed 
line), 17.0 (dotted-dashed line), 20.0 (dotted line) and 22.5 (full 



line). The dotted vertical lines, passing by L 



(2) 



0.78 m~^ and 



.(2) 



"^, are plotted as a visual guide. 



FIG. 1: Plot of S2 = 1121 (I/) as a function of mL. 
Recall that, in the present case, Eq. ( fT9] l becomes 



(2)/ 



92i{L,X) 



X 



i + An2i(i)' 



(23) 



The divergence of 1121 (L) as i —s^ ensures that, inde- 
pendently of the value of A, g2i{L, A) approaches in this 
limit and, therefore, the system presents a kind of asymptotic- 
freedom behavior for short distances. On the other hand, 
since 1121 (£) assumes negative values for L > i,„ijj, the de- 



nominator of Eq. (I23l l will vanishe at a finite value of L if 
A is sufficiently high. This means that, starting from a low 
value of L (within the region of asymptotic freedom) and in- 
creasing the size of the system, 921 will diverge at a finite 
value of L, Lc (A), if A is greater than the "critical value" 

A^^^ = (-n™")-i ~ 22.47. We interpret this resuh by Stat- 

(2) 
ing that, in the strong-coupling regime (A > Ac ) the system 

(2) 
gets spatially confined in a segment of length Lc (A). The 

behavior of the L-dependent coupling constant as a function 

of mL is illustrated in Fig. [21 for some values of the fixed 

coupling constant A. 

It should be emphasized that we are treating the massive 
GN model with an arbitrary (but fixed) fermion mass. In this 
case, the model does not possesses chiral symmetry which is 
explicitly broken. Since we do not expect that this symme- 
try appears beyond the critical value (with the radiatively cor- 
rected mass vanishing identically), the instability indicated by 
the divergence of the coupling constant is interpreted as sig- 
naling a spatial confining transition. 

For A = Ac , by definition, the divergence of g2i{L, A) is 
reached as L approaches the value that makes 1121 minimal. 



-(2) 



In the other hmit, since 



which we have denoted by L„ 

921 {L,\ ^ cx)) = n2U 

of 1121 (L), as A ^ 00. In other words, the confining length 



l(L), ii^^(A) tends to l[^1^, the zero 



SL, when A 



Lc"' (A) decreases from the maximum value L 

(2) (2) 

Ac , tending to the lower bound Li,/„ in the limit A 

^ ' ^ min 

The behavior of Lc 
later. 



as a function of A, will be presented 



B. Compactified 3-D GN model at T = 

For the 3-D model at zero temperature with two com- 
pactified dimensions (d — 2), denoting the compactification 
lengths associated with the two spatial coordinates xi and X2 
by Li and L2 {mT^ j \fb\ and mT^ j \fb2, respectively), for- 
mulas in Eq. (fTsT i and Eq. ( fTSl ) give 



1132(^1, &2) 



^/27r^ 



[2L^3i(2;6i)-L^3i(2;4fei) 



+ 2F3i(2; 62) - F3i(2; 4&2) - 2F3i(l; 62) 

- 2F3i(l; 61) + F3i(l; 4&i) + F3i(l; 462) 
+ 4F32(2;6i,62)- 2^^32(2; 461,62) 

- 2F32(2; 61, 462) +L^32(2; 461,462) 
-4F32(1; 61, 62) + 2i^32(l; 461,62) 

+ 2^32(1; 61, 462) - F32(l; 461,462)] .(24) 

The functions F-^j (j = 1, 2), specified in Eq. ( fT4] l. involve the 
Bessel functions of order ±i, which are expressed in terms of 
elementary functions: 



K^Az) = ^ 



exp(— z) 



2z 



(25) 



Thus, the series defining the functions F3J, for both j — 1,2, 
are geometric series which can be summed up. Using Eq. (l25t 
in the expression of Eq. (fT4] i and replacing 6; by L^"^ (which 
corresponds to taking all the compactification lengths mea- 
sured in units of mT^ , as will be done from now on), we ob- 



tain 



1132(^1,^2) — 7^ 



m 

2^ 



— log(l + e-^i; 
^1 



1 



1 + 6-^1 
1 



+ — log(l + e ^■') + - 
L12 i 

in 

+ -[G2{Li,L2)-2G2{Li,2L2) 

-2G2(2Li,i2) + 4G2(2Li,2i2)] 

where the function G2 {x, y) is defined by 

00 
G2{x,y) = y^ cxp (-y^x^n^ + y^f 

1 



(26) 



ra,i=l 



^ x^v? + y^Z^ 



(27) 



The numerical computation of 1132 (£1, ^2) is greatly facili- 
tated by the fact that the double series defining the function 
G2(y, 2:) is rapidly convergent. The need of truncating the 
summations at a larger value n — M when L^ and L2 are 
very small, as in the case of 13 = 2, still exists but to a much 
less extent. 

It is to be noticed that if either compactification length L\ 
or L2 tends to 00, all terms depending on it disappears from 
Eq. ( |26] l and we regain the finite bubble diagram for the case 
where only one spatial dimension is compactified in the 3-D 
model [35] . Now, if both L\ and L2 tend simultaneously to 
oo^ 1132 goes to zero and 1732 -^ A, as expected. On the other 
hand, if either L\ or L2 tends to 0, 1132 -^ +00 implying that 
the system gets asymptotically free, with the effective cou- 
pling constant vanishing in this limit. The overall behavior 
of the bubble diagram is illustrated in Fig. |3] where we draw 
the contourplots of n32(£i, L2)/m. The full line in Fig.[3]is 
the locus of the points such that 1132(1/1, L2)/m = 0, which 
for large Li (L2) approaches the straight line Li — l.lAm'^ 
(L2 — 1.14 771^^); n32(£i,L2) is positive below this curve, 
negative above it, and reaches an absolute minimum, 1132" — 
-0.00986 m, at the point Li = L2 ~ 2.10 m-\ 

The fact that 1132 assumes negative values in the whole re- 
gion of the parameter space {Li,L2) above the full line in 
Fig. [3] implies that, for large enough values of A, (732 will di- 
verge at finite values of Li, i = 1,2. However, to avoid un- 
necessary complication, our analysis is restricted to the case 
where the system is confined within a square of size L, by con- 
sidering Li = L2 = L. In other words, we shall concentrate 
on the behavior of 1132 along the diagonal of Fig. [3] We plot, 
in Fig.m 5*3 (L) = 1132 (i)/m as a function of i. We find that 
S{L) vanishes for a specific value of L, L[J^, being negative 
for all L > i„jj,j. Also, 1132 (L) reaches an absolute minimum 

(3) 

(negative) value for a value of L we denote by Lmax- We find, 
numerically, that L^^l^ — 1.30 m~^ and Lmax — 2.10 m~^ , 
with 1132" — — 0.00986 r?T,, as stated before. This behavior of 
1132 (i) has profound implications on the effective coupling 
constant. 



\ \ 



1.5 2 2.5 



3.5 4 4.5 



FIG. 3: Contourplots of 1132(1/1, 1/2)/f7i, with Li and I/2 in units of 
m^^. The open dashed line corresponds to 1132(^1, L2)/m = 0.04, 
the full line gives the points where 1132 — 0, while the closed curves 
are for negative values of n32/m, -0.0091, -0.0095 and -0.0098 
(dashed-dotted, dashed and dotted lines, respectively). The dot is the 
location of the absolute minimum of 1132(1/1, 1/2), which occurs for 
Li= L2~ 2.1 m"\ 




FIG. 4: Plot of 53 = 1132 (L) /m as a function of L, in units of 7 



With D = 3 and d = 2, Eq. (O is rewritten as 

532 (i, A) 



1 + An32(i)' 



(28) 



i(3) 



'^'"Y 



101.42 m- 



and we find that, for A > Ac 

the denominator in Eq. (l28l l will vanish for a finite value 
of L, Lc (A), leading to a divergence in the effective cou- 
pling constant. The behavior of the effective coupling con- 
stant as a function of L, for increasing values of the fixed 
coupling constant A, can be illustrated showing the same 
pattern as that of Fig. |2] for the preceding case. Similarly, 
we find that the divergence occurs at Lc (A) which satisfies 

Again, we interpret such a result 



(3) 



(3)/ 



< Lc (A) < Lmax- 



(3) 
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by considering the system spatially confined in the sense that, 
starting with small L (in the region of asymptotic freedom), 

the size of the square cannot increase above Lc (A), since 

c/32(L, A) ^ oo as L ^ Lf\\). 



C. D — 4, case at zero temperature 



In the case of the A-D GN model with all three spatial co- 
ordinates compactified, replacing bi by L^^ (again, Li mea- 
sured in units of m^^) into Eqs. ( fT3] l. (O and ( fTSI l gives 

n43({ij) = m^2Hi{2Li) + 2Hii2L2) + 2Hi{2L3) 

- Hi{Li) - Hi{L2) - HiiLs) 
+ 2H2iLi,L2) + 2H2iLi,L3) 
+ 2H2{L2,L3)-AH2iLi,2L2) 
-AH2iLi,2L3)-AH2i2L,,L2) 
~4H2{2Li,L3)-AH2{L2,2L3) 

- AH2{2L2, L3) + 8H2{2Li,2L2) 
+ 8^2(2X1, 2L3) + 8H2{2L2, 2L3) 



-4H3{Li,L2,L3) + 8H3{2Li,L2, 
+ 8H3{Li,2L2,L^) + 8H3{Li,L2, 
-16i/3(2ii, 2^2,^3) 
-16i/3(2ii,i2,2L3) 
-16ff3(ii,2L2,2L3) 
+ 32^3(2^1, 2L2,2L3)], 

where the functions Hj, j = 1,2,3, are defined by 

Ki{xn) 



2L3) 



Hi{x) 
H2{x,y) 



2tt 



TT^ ^ — ' ixn) 

n—l '- ^ ' 

_, 00 

-2 Y. [^0 (y^^^^^TW^) 



(29) 



(30) 



,,i=i 



Kx{\Jx'^n^ + y'^P) 



H3ix,y,z) 



1 

2^ 






(31) 



E 

i,l,r=l 



+ z^r' 



Kiiy/x^n'^ + yn^) 



(32) 



Results for one or two compactified dimensions are ob- 
tained from Eq. ( |29] | if two or one of the compactification 
lengths become infinite. For example, taking ^2,^3 —> 00, 
we get n4i(Li) — m}\21I\(2L\) — iI\(L\)\, which assumes 



the minimum value IIJJ^'" 

2.01 niT^, leading to the critical value A 



-0.001866 m^ at L = L 



(41) 



(41) 



535.91m- 



and vanishes at L 



(41) 



1.43 771 ^. The analysis made 



in Ref. fSS"] can be extended to this case of the 4-dimensional 
space with only one spatial coordinate compactified. Simi- 
larly, we could discuss the situation with two compactified 
spatial dimensions. 



Here, instead of dealing with all possibilities, we concen- 
trate on the case where the system is confined to a cubic box, 
that is, we take Li ^ L2 — L3 = L. With equal compactifi- 
cation lengths, Eq. ( |29] l becomes 

n43(i) = m'^[6Hi{2L)^SHi{L) + 6H2{L,L) 

- 24:H2{L, 2L) + 24iJ2(2L, 2L) 

- AH-i{L, L, L) + 24iy3(L, L, 2L) 

- 48H3{L, 2L, 2L) + 32i?3(2L, 2L, 2L)] . 

(33) 

This quantity is plotted in Fig. |5] which shows that it has 
the same behavior as its counterparts for D = 2 and D = 
3. We find numerically that 1143 (L) vanishes for L = 



(4) 



1.68 771 ^, being negative for L > L^l^, and as 



L 

sumes the minimum value, 1143' 



-0.0022751 m2, when 



-(4) 



0.02 



0.015 



0.01 



0.005 



-0.005 



2.37TO- 



FIG. 5: Plot of 54 = Il43{L)/m'^ as a function of L, in units of 

»v,-l 



As in the other cases, the large- A^ coupling constant. 



543 (^, A) 



A 



i + An43(i) 



(34) 



diverges at a finite value of L, Lc (A), if A > Ac — 
— (1143")^^ ~ 439.54?7i~^, meaning that the system gets 

confined in a cubic box of edge Lc (A) which is bounded 

in the interval between the values i,^;,! and L,nax- 



D. Dependence of Li ' on A 

For the cases we have analyzed above, namely D = 2,3,4 
with all spatial coordinates compactified and Li ~ ■ ■ ■ = 
Ljj^i — L, we find that the confining length Lc (A) lies 
in a finite interval, 



iPnA)e(4^^,L(f.L 



(35) 
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where the maximum value corresponds to Ac , the minimum 
value of the fixed coupling constant A allowing spatial con- 

-(D) 

question emerges of how Lc (A) changes as A increases from 
Ac to infinity. 

For a given value of A (> Ac ), the confining length 
Lc (A) can be found numerically by determining the small- 
est root of the equation 



feD-i(i,A) = -[l 



An 



DD-l 



(L)]=0. 



(36) 



That is, following the interpretation provided before, start- 
ing from small values of L, the first value at which g]j]j_i 
vanishes does provide the confining length of the system. 



Lc (A). In Fig. |6l we plot Lc (A) as a function of I = 
A/aP', for the cases D = 2, 3, 4. 




10 12 



FIG. 6: Plot of the confining length (in units of m ^), as a function 



AD) 



of I = X/K ', for D = 2, 3, 4 (dashed, dashed-dotted and full lines 
respectively); the horizontal dotted lines correspond to the limiting 



.(D) 



-(D) 



values I/Jjji„ and Lmax (given in the text), plotted as a visual guide. 



IV. COMPACTIFIED GN MODEL AT FINITE 
TEMPERATURE 

We shall now consider the effect of raising the temperature 
on the effective coupling constant for the GN model with all 
spatial dimensions compactified. Finite-temperature effects 
are introduced through the compactification of the time co- 
ordinate, with the compactification "length" given by L/j = 
/3 = l/T, where T is the temperature. It is important to em- 
phasize that, although in an Euclidean theory time and space 
coordinates are treated on the same footing, the interpretation 
of their compactifications are rather distinct; while compact- 
ification of spatial dimensions can be thought as describing 
confined fields, time compactification corresponds to taking 
the system in thermal equilibrium at temperature /3^^. 

We generally expect that the dependence of II^i^i and g^D 
on /3 should follow similar patterns as that for the dependence 
on L. Like in the case where any compactification length Li 



tends to zero, we find that lioDiL^ /3) — )■ cx3 as /3 — >■ (T — >■ 
oo), implying that gDD -^ independently of the value of 
the fixed coupling constant A. This means that we have an 
asymptotic-freedom behavior for very high temperatures. 

For /3 — > oo (T — > 0), the behavior of IiDD-i{L) has been 
described earlier: for sufficiently high values of A, the sys- 
tem is confined in a (13 — 1) -dimensional cube of edge Lc ■ 
Based on these observations, we expect that, starting from the 
compactified model at T = with A > Ac , raising the 
temperature will lead to the suppression of the divergence of 
goD and the consequent spatial deconfinement of the system, 
at a specific value of the temperature, T^ . The way such a 
"deconfining" transition occurs and the determination of the 
deconfining temperature for the cases of 13 — 2, 3, 4 are the 
points addressed in the next subsections. 



A. 2-D compactified GN model at T / 

To account for the effect of finite temperature on the 2-D 
compactified GN model 115911 . we take the second Euclidean 
coordinate (the imaginary time, X2) compactified in a length 
L2 — P = l/T. In this case, replacing 61 = L^'^ and 
62 — /3~^ {L and /? measured in units of m"^) into Eqs. ( fT3] ), 
(fT4l) and ( fTSI ). the L and /3-dependent bubble diagram can be 
written as 

n22(L,/3) - 2E^{2L)-E^{L) + 2E^{2p)~E^{li) 
+ 2E2iL,(3)-4E2{2L,(3) 
-4S2(i,2/3) + 8^2(2L,2/3), (37) 

where the function Ei (x) is given by Eq. ( |22] | and the function 
E2{x, y) is defined by 



ri,/=l 

+ (^/x'^rfi +y'^lA Ki (^/x'^'n? + y'^P 



(38) 



We first notice that, due to the behavior of the Bessel func- 
tions Ko{z) andKi{z), all /3-dependent terms in Eq. (|37]i van- 
ish in the limit /3 — > cx) and so 1122 (i, /?) reduces to the ex- 
pression for zero temperature, 1121 (L). On the other hand, if 
/3 — > 0, 1122 {L, P) —^ 00 and, independently of the value of A, 
the system becomes asymptotically free. Therefore, we expect 
that, within the strong-coupling regime, raising the tempera- 
ture leads to the suppression of the divergence of 322 and the 
disappearance of the spatial confinement. In other words, for 
a given value of A > Ac , there exists a temperature, T^"' (A), 
above which (722 has no divergence and the system is spatially 

deconfined. 

(2) 
The deconfining temperature T^ (A) is determined by an- 
alyzing the behavior of 17^2^ (L, (3, A) as T is increased. This 
process is illustrated in Fig. |7] where we plot 17^2^ (L, /3, A) as 
a function of L, for some values of (3 and a fixed value of 
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A > Ac . For this example with A = 30, we find that the min- 
imum value of g22 vanishes for [i = /3^ ~ l.lSrn"^ and is 
positive for P > P^^ . Thus, the deconfining temperature, for 

A = 30, is given by rf ' = Wd^y^ - 0.87 m. The full de- 
pendence of deconfining temperature on A will be discussed 
later. 




FIG. 7: Inverse of the effective coupling constant g^^, with A = 30 
fixed, as a function of L (in units of m~ ), for some values of /3 
(in units of m~^): 2.4, 1.15 and 1.0 (dashed, full and dotted lines, 
respectively). 



where G2{x,y) is given by Eq. dZTJ ) and the function 
G3{x,y, z) is defined by 



G 



3{x,y,z) — 2. exp I — Y'a^^fi^ + y^Z^ + 2^r2 



,l,r=l 



^ x^v? + j/2p + z^r^ 



(40) 



Notice that, taking fi -^ co, Eq. ( [39] l reduces to 1132 (i), ob- 
tained from Eq. ( |26] | with Li = L2 = L. 

As before, the increase of the temperature destroys the 
spatial confinement that exists for A > Ac at T = 0. 
We can determine the deconfining temperature by search- 
ing for the value of /3(A) for which the minimum of the 
inverse of the effective coupling constant, g^^{L,(3,X) — 
(1 + An33(L, 13))/ X, vanishes. For example, taking the spe- 
cific case of A = 110 ?ti^^, we find /3^ ~ 1.65 m^^ which 



.(3) 



corresponds to the deconfining temperature of T^ ~ 0.61 m; 
this result can be illustrated in a figure with the same pattern 
as that appearing in Fig.|7]for the D = 2 case. 



C. D = 4 case at T / 



B. 3-D compactified GN model at finite T 

We now take the time coordinate {X3) compactified in a 
length /? = 1/T to investigate the temperature effect in the 
3-D compactified GN model. In Ref. f3o] the 3-D model, 
with only one spatial dimension compactified, was treated at 
finite temperature 1I6OII . Here, we deal with the fully compact- 
ified model. Taking 61 = 62 = L^^ and fixing 63 — /3^^ in 
Eqs. ( fT3] l, (O and ( fTSl ), the L-/3 dependent bubble diagram is 
given by 



n33(^,/3) 



m 

2^ 



L 



log(l 



log(l + e-^) 
[G2{L,L) + 



' l+i 


,i 


1 




l + e/3 




2G2(i,/3) 


-AG2{L,2L) 



4G2(2i, /3) - 4G2(L, 2/3) + 4G2(2L, 2L) 

-8G2(2i,2/3)-2G3(i,£,/3) 

-4G3(£,i,2/3)+8G3(2i,i,/3) 
8G3(2L,2L,/3)-16G3(2L,L,2/3) 
16G3(2L,2i,2/3)], (39) 



ForthefuUy compactified 4-£'GN model, fixing &i = &2 = 
63 = L~^ and 64 = /3~^ {L and /3 in units of m^^), we find 
from Eqs. O, (O and ^ that 

n44(L,/3) = m'^[QHi{2L)-?.Hi{L)+2Hi{2f3) 

- Hi{(3) + 6H2{L, L) + 6H2{L, /3) 
-24i/2(i,2L)-12iJ2(i,2/3) 

- 12H2{2L, /3) + 24H2(2L, 2L) 
+ 24H2{2L,2(3) - m3{L,L,L) 

- UHsiL, L, /3) + 24:H3{L, L, 2L) 
+ 48H3{L, 2L, /3) + 24iJ3(i, L, 2/3) 

- 48iJ3(i, 2i, 2L) - 48H3{2L, 2L, /3) 

- 96H3{L, 2L, 2/3) + 32iJ3(2i, 2i, 2i) 
+ 967l3(2i, 2L, 2/3) + 8i/4(i, i, i, ^) 

- 48F4(L, L, 2L, (3) + 192Hi{L, 2L, 2L, (3) 

- l&Hi{L, L, L, 2/3) + 96H4,iL, L, 2L, 2(3) 

- MHi{2L, 2L, 2L, /3) 
-192Hi{L,2L,2L,2(3) 

+ 128Hi{2L,2L,2L,2(3)], (41) 

where the functions Hi, H2 and H3 are given by Eqs. (|30II32| |, 

and H4{x, y, z, w) is defined by 
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^ oo 



■ ,/,r,s— 1 



Kq I \J x^v? + ip-X^ + z^r^ + w^s^ ) 






(42) 



r 



Proceeding as before, we determine the deconfining tem- 
perature by searching for the value of /9(A) for which the 
minimum of the inverse of the effective coupling constant, 

g44^(L,/3,A) = (1 + An44(L,/3))/A, vanishes. For ex- 
ample, taking the specific case of A = 620 m^^, we find 

fi\ ~ 1.707™^^ which corresponds to the deconfining tem- 



perature of T^ 



(4) 



0.59 m. 



(D) 



D. Dependence of T\ on A 

We now determine the dependence of the deconfining tem- 
perature on the fixed coupling constant. For the system with 
all spatial dimensions compactified and at finite temperature, 
we have 



g^}3(L,/3;A) = i[l + An 



DD 



(i,/3)] 



(43) 



As discussed before, the system is deconfined at a given tem- 
perature if the minimal value of g^^(L, /3; A) with respect to 
changes in L is positive. It was also argued that, no matter 
how high the value of A is, the system becomes deconfined 
above a given temperature, T\ (A). In fact, one expects that 



^d y'^l ^ -^min' -^ ma> 



(44) 



SD) 



with the limiting values corresponding to Ac and A — > 
respectively. 

From Eq. ( |43] l, we find that minj^} ^^^(L,/?; A) 
[1 + A Afi3(/3)] /A, where we have defined the function 

Afi5(/3)=minni3i3(i,/3). 



oo. 



(45) 



For a fixed value of A, the behavior of the minimum value 
of g]j]j{L,(3; A) relative to changes in L is dictated by the 
function Md{(3). Notice that, Md{13 ^ cx)) = ng'™^^ = 

-[Ai^V^ while Md{P) ^ oo as /3 ^ 0. 

Consider, initially, the case D — 2. The function Af2(/3) 
is illustrated in Fig. |8] We find that, as /3 decreases from oo 
(i.e., T increases from 0), the minimum value of 1122 (i, /3) 
(with respect to changes in L) starts to decrease from nega- 
tive values, passes through the lowest value and then starts to 
increase, reaching zero at a certain value of j3 below which 

('2) 

the minima of 1122 are positive. Thus, for A = Ac , increas- 
ing the temperature from zero, the system remains confined 



until the temperature reaches the value T^ 



(2) 



[/3, 



(2) 



0.65 m, corresponding to the finite solution of the equation 



1 + xf^NhiP) = 0, /3 



g(2) 
Mniax 



1.54 m ^, which is in- 



dicated by the vertical dotted line in Fig. [8] Now, if we take 



A > Ac , the deconfining temperature is obtained from the 
solution of the equation 1 + AM2 (/3) = which can be deter- 
mined from the intercept of the horizontal line at — A ^ ^ (lying 
above the line —[Ac ]~^ and below the j3 axis) and the graph 
of A/2(/3). Naturally, as one takes A — ?► 00, the existence of 
a solution of this equation requires Af2 (/3) — > 0; this lower 



(open) limit occurs at the value /? 
responding to the temperature T, 



(2) 



"nun 



0.776 m ^, cor- 



t 1.29 m. The depen- 

.(2), 



dence of the deconfining temperature Tj "^ on A is determined 
numerically. We find that, for A not close to Ac (that is, for 
A > 1.5 aI^^), ^^^^(A) = L^c\\) within six decimal places; 
this approximate equality, valid for large values of A, is a con- 
sequence of the symmetry of the expression for 1122 {L, (3) by 

/ON 

the change L -^ j3. The whole behavior of Tj ^(A) is shown 
in Fig.[TOl together with the other cases. 



Klili) 



FIG. 8: Minimal values of 1122 (i; /3) with respect of changes in L, 
as a function of /3 (in units of m~^). The dotted horizontal line 
corresponds to the value —[Ac ]~^- 



It should be remarked that the absolute minimum of 
n22(L,/?), l^r - -0.0466, is sHghdy smaller than H^Y" 
and occurs at the point with L == /? ~ 1.98 m^^. This means, 
as shown in Fig.[8l that the minimum value of AhiJi) does 
not occur at zero temperature. This leads to an anomalous sit- 
uation, for a small range of values of the coupling constant A 
(21.46 < A < X^c^ ~ 22.47), in which no singularity exists 
at T = but the equation 1 + AA/2 (/?) = possesses two so- 
lutions for finite values of /?; this would imply that the system 
being unconfined at T — would get confined at a finite tem- 
perature and then becomes unconfined again at a smaller value 
of /3. Such a situation, which emerges from the mathematical 
structure of the zeta-function regularization, has no physical 
meaning and will be discarded; we shall only consider the 
strong coupling regime which is free from pathologies. 

For D — 2>, the graph of the function Mj,{(i)/m has 
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a form similar to that of Af2(/3) (illustrated in Fig. [8]l so 
that the same reasoning leads to the deconfining temperature 

T^^\X) = [/3'f\x)]-\ where;3^^^(A) is the finite root of the 
equation 1 + AM3(/3) = 0. In this case, the limiting val- 
ues are /3max — 1.85 rn^^ and (3j^(^ — Llim^^, associated 

(3) (3) 

with A = Ac and A — > oo, corresponding to T^^^^ — 0.54m 
and T^L ^ 0.88m. The plot of T^^\\), as a function of 

I — A/ Ac , is also shown in Fig. [TO] Distinctly from the 
D = 2 case, where all the expressions are symmetric by the 
change L o /3, here we find that ^^^^A) < L^^^(A) for all 
A > Ac , the difference being of the order of 10%. The maxi- 

f3) 

mum value of T\ (~ 0.87 m), which occurs as A ^ cx), cor- 
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responds to /3^^^ reaching the value of the minimal confining FIG. 9: Minimal values of 044 (i; l3)/m^ with respect of changes in 



length (~ 1.14 m ^) when only one coordinate is compacti- 
fied. 



L, as a function of /? (in units of m ^). The dotted horizontal hues 
correspond to the values —[Ac ]~^ < Afmin < Mmax. 



The D = 4 case is more subtle due to the behavior 
of the function M4{fi), which is presented in Fig. |9] We 
find that the minimum value of 1144 (L,/3) (with respect to 
variations of L) starts to increase from negative values as 
/3 is diminished from cx), reaches a local maximum value 
(-^^max — — 0.0018356 m^), decreases to a local minimum 
(Afniin — —0.002091 m^) before increasing to reach zero and 
become positive. Therefore, for A = Ac = — [A/4(/3 — >■ 
cx))]^^ ~ 439.54 m~^, no finite solution of the equation 

1 + Ac M4(/3) — exists, which means that the system is 
deconfined if the temperature is greater than zero, no matter 

(4:) 

how small it is; that is, for A = Ac , spatial confinement is 
only possible stricdy at T = 0. For Ac''' < A < 478.2 m~^ 



(- 



-il/jjjjjj), the equation 1 + XM^^P) = possesses one 



solution occurring at finite /3, which could eventually be in- 
terpreted as leading to a deconfining temperature. However, 
if 478.2 m"2 < A < 544.8 m"^ (^ -A/^-i^), the equation 
1 + XM4{(3) — has three distinct finite solutions. If we in- 
terpret the highest one as giving the deconfining temperature, 
we would have to face the puzzling situation in which the sys- 
tem would reenter a spatially confined phase for /3 ranging 
between the other two smaller solutions. Such an anomalous 
behavior can be avoided if we redefine the strong coupling 
regime by considering the range —M~l^ < A < 00. 

With such a redefinition of the strong coupling regime for 
D — 4, which amounts to considering the lowest value of A 



leading to spatial confinement as being A 



(4) _ 



-M- 



544. 



we find a well defined deconfining temperature 



obtained from the intercept of the horizontal line at — A^^ and 
the curve Af4(/3). We get (3^L ~ 1.80 m~^ (indicated by the 

vertical dotted line in Fig. |9]l and /?,„;„ — 1.43 m^^ (where 
M4 vanishes), corresponding to the limits of the deconfining 

0.70 m., respectively. 



temperature T^-^^ 



0.55 m and Tifax 



(4) 



Similar to the D = 3 case, the value of /3,„4, giving the upper 

bound for Tj (A — > 00), is identical to the smallest confining 
length when only one spatial dimension is compactified. The 
overall behavior of TJ (A), found numerically, is presented 
in Fig. [TO]together with the cases of £> = 2 and 3. 
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FIG. 10: Deconfining temperature T^ (A) (in units of m), as a 
function of / — A/ Ac , for D — 2,3,4 (dashed, dashed-dotted 
and full lines respectively); the horizontal dotted lines correspond to 
the limiting values T^^^ and Tmax (given in the text), plotted as a 
visual guide. 



V. CONCLUDING REMARKS 

We have used methods of quantum field theory in toroidal 
topologies to investigate the behavior of the D-dimensional 
massive, A^-component, Gross-Neveu model with compact- 
ified spatial dimensions. The model is treated both at zero 
and finite temperatures. We calculate the large- A^ coupling 
constant, go, as a function of the compactification length L, 
the temperature T = /S^^ and the fixed coupling constant A, 
i.e. gjj = goiL, /3;X). We find that, for either L ^- or 
T -> 00, irrespective of the value of A, go tends to indicat- 
ing that the system presents a sort of asymptotic-freedom be- 
havior in these limits, where the effective interaction between 
the fermions vanishes. 

For T = 0, in the strong-coupling regime (A > Ac ), in- 
creasing L from low values (within the asymptotic freedom 
region) leads to a divergence of gjj at a finite critical value 
L — Lc (A), signaling the existence of a second-order phase 
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transition, as suggested by the linear response theory. Also, 
since we consider a four-fermion interaction model at zero 
chemical potential, we should not expect a first-order phase 
transition for any value of the parameters characterizing the 
system, as happens for massless models. We interpret this sin- 
gularity as indicating that the system gets spatially confined in 
a (Z? — l)-dimensional box of edge Lc (A). We have shown 
that Lc (A) = foWrn^^, where the functions /13(A) are 
plotted in Fig.|6l for D = 2, 3, 4. 

As T is raised from 0, with A > Ai , the minimum value 
of g]j^{L, (3; A), with respect to changes in L for fixed val- 
ues of /3, increases from negative values reaching at a tem- 
perature T^ (A) above which g^, does not present any di- 
vergence. We interpret this fact as the system being decon- 
fined for temperatures higher than T^ (A). To avoid any 

anomalous behavior, we have redefined the value of Ac aug- 
menting the lower bound which defines the strong-coupling 
regime when D ~ 4,. In any case, we find that TJ '(A) = 
h£){X) m, where h]j{X) are the functions plotted in Fig. (TO] 
for D — 2,3, 4. It is to be noted that, for the redefined value 
Ac ~ 544.8 m^^, the zero-temperature maximum confining 



length is given by L 



(4) 



2.00 m" 



(D) 



It is worth emphasizing that the dependencies of Lc and 
T^ on the parameters A and m are intrinsic results of the 
model, that is, they do not emerge from any adjustment. The 
dependence on m is precisely what one expects from di- 
mensional arguments, with L and /? being proportional to 
m^^ in natural units. But since the functions /d(A) and 
/i£)(A) take on values in finite intervals (in all cases, con- 
tained in [0.5,2.5]), we find that extremely light fermions 
(to -^ 0) are not confined at all, while extremely heavy ones 
(m —}' 00) would be strictly confined in a dot, no matter what 
the value of A. Also, it should be noticed that the product 
Vd{X) = Li^\x)T^'^\x) = /d(A)/id(A) is dimension- 
less and very close to the unit in the whole strong coupling 
regime, A > Ac , irrespective of the value of D; actually, 
one finds that Vd{X) G (1.00, 1.18) for all cases. Such an 
inverse relation between the confining length and the decon- 
fining temperature is what one would expect from strong in- 
teraction and QCD physics; if the length of confinement of 
a fermion is small, the energy required to overcome its con- 
fining barrier is large, and vice-versa. Roughly speaking, one 



has a sort of uncertainty relation L ~ 1/p, with p being the 
momentum which, for relativistic particles, is proportional to 
the energy. Anyhow, the questions of how the fermions are 
confined and how they get unconfined are not answered with 
our analysis. 

Since we have completely determined the relevant depen- 
dence of Lc and T^ on A, for any value of A (> Ac ), 
estimates of values of the confining length and the deconfin- 
ing temperature require the specification of the mass of the 
fermions. Viewing the GN model as an effective model for 
strong interaction, a natural choice is to take the constituent 
quark mass, m ss 350 MeV ~ 1.75 fm~^ 1I6II1 . With such a 
choice, we obtain the limiting values of the confining length 
and the deconfining temperature presented in Table 1 . 



TABLE I: Limiting values of L^^' (A) [fm] andrj^'(A) [MeV], for 
m fa 350 MeV ~ 1.75 fm"\ 



D 


r(D) 




min 


-l- max 


2 


0.45 


0.96 


227 


451 


3 


0.74 


1.2 


189 


304 


4 


0.96 


1.14" 


193" 


245 



"These values of Ljnax and T ■' con'espond to the redefined value of AJ, . 

We find that the range of variation of Lc and T^ % as 

A runs in the strong coupling regime (> Ac ), is relatively 
small and decreases as D increases. These values compare 
amazingly well with the size of hadrons (e.g., the experimen- 
tally measured proton charge diameter is w 1.74 fm [62J) 
and the estimated deconfining temperature (« 200 MeV) for 
hadronic matter [631. 

To summarize, we have established the existence of a 
phase transition in the massive large-A^ GN model with 
compactified spatial dimensions. It remains to be investigated 
other aspects of this transition like effects of finite density 
and critical exponents. Such a study is left for a future work. 
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